Mechanism of vortices generation for a trapped superfluid under oscillatory excitation 
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The vortex nucleation and the emergence of quantum turbulence induced by oscillating magnetic 
fields, introduced by 0,01, left a few open questions concerning the basic mechanisms causing those 
interesting phenomema. Here, we report the experimental observation of the slosh dynamics of a 
magnetically trapped *^Rb Bose-Einstein condensate (BEC) under the influence of a time-varying 
magnetic field. We observed a clear relative displacement in between the condensed and the thermal 
fraction center-of-mass. We have identified this relative counter move as an out-of-phase oscillation 
mode, which is able to produce ripples on the condensed/thermal fractions interface. The out-of- 
phase mode provides good evidences of a mechanism involved in the vortex nucleation and evolution 
when excited by time dependent magnetic fields. 

PACS numbers: 47.37.-(-q, 67.10. Jn, 67.85.-d 
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I. INTRODUCTION 

The research evolution involving atomic superfluid 
during the last decade has been quite significant. It 
started with the demonstration of the superfluidity as- 
sociated with Bose-Einstein condensation (BEC) 3^, the 
nucleation of vortices , the formation of vortex lattices 
0, and finally the first observation of quantum turbu- 
lence in a ^^Rb BEC 0- Such evolution presents the 
full potential of using atomic superfluid samples as work- 
benches for many phenomena that could only be exper- 
imentally accessed before using liquid helium Q. Re- 
cently, an oblate sample of trapped atoms, corresponding 
to a 2D-superfluid, has been explored in many aspects 
related to superfluidity. A repulsive obstacle passing 
through this 2D system has generated vortex-antivortex 
pairs and its rotation dynamics could be studied [7]. In 
the meanwhile, theory has also progressed in many di- 
rections [8l-[l0|. In a recent experiment on generation of 
vortex by applying an oscillatory field in a trapped Bose- 
Einstein condensate [1] , the proliferation of vortices with 
amplitude and time of excitation [iTl [l^ results in a con- 
figuration of vortex filaments in many spatial directions, 
allows to obtain a turbulent cloud of atomic superfiuid 
Ql- In those already reported experiments, the mech- 
anisms from which vortices were generated was still an 
open question. Ex plan ation were based on a few the- 
oretical simulation [13|. Experimental evidences of nu- 
cleation of vortices in oscillatory conditions, previously 
to vortices proliferation and turbulence, was a missing 
fact in the sequence of experiments. This report brings 
the first evidences of excitations that can promote vortex 
nucleation during oscillatory excitation. 

To understand our present experiment we shall start 
with a widely explored scenario in liquid helium, below 
the A-point, which corresponds to the two-fiuid system 
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created by the coexistence of a superfluid and a nor- 
mal fluid fraction In this case, the normal fluid 
can be regarded as a classical viscous fluid. Experiments 
in which a counterflow between the two parts became 
possible, resulted in many important experimental inves- 
tigations including the equivalence of classical flow using 
the Reynolds number criterion [l^ and the generation of 
vortices. 

The two-fluid system seems appropriate to help an- 
swering questions concerning the onset of dissipation on 
the superfluid flow. One of most interesting counterflow 
experiments is the thermal counterflow . In that case, 
the temperature gradient generates a superfluid flow op- 
posite to the normal fluid flow, giving rise to a mutual 
frictional force between the two fluids [16]. Although 
thermal counterflow was not yet been realized in a sample 
of trapped atomic superfluid, a collective mode present- 
ing an out-of-phase oscillation between the condensate 
and thermal cloud has been investigated systematically 
by [13,0. 

In this report we show that, under the influence of a 
oscillatory magnetic field, a trapped BEC may evolve to 
a scenario where the condensed and the thermal frac- 
tions present a relative out-of-phase motion, character- 
izing the onset of counterflow (second sound type exci- 
tation). Such situation may be useful to investigate the 
circumstances in which quantized vortices and turbulence 
take place during oscillatory excitation. 



II. EXPERIMENTAL SETUP AND TIME 
SEQUENCE 

The experimental sequence used to produce the BEC 
and nucleate vortices in our condensate can be summa- 
rized as follows. First, a ^^Rb BEC containing about 
1 X 10^ atoms with a small thermal fraction is produced in 
a cigar-shaped QUIC magnetic trap with: i/r = 210 Hz, 
and Vz — 23 Hz trapping frequencies. While the BEC 
sample is realized and still held inside the trap, a weak. 
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sinusoidal, time- varying magnetic field is turned on. This 
field is superposed to the QUIC trap and is generated by 
an extra pair of Helmholtz coils mounted with the sym- 
metry axis aligned in a small angle relative to the lofFe 
coil axis. Since the coils are axially misaligned, the gra- 
dient near the trapping region present field components 
parallel to each of the trap eigen axis. This time-varying 
approach was originally developed to investigate coherent 
mode excitations in our BEC samples. 

A combination of shape shifting, trap minimum dis- 
placement and axial twisting move couples to the atoms. 
The perturbing field was on for time lapses usually rang- 
ing from 20 to 100 ms; and it oscillates sinusoidally in 
time with frequencies in between 150 to 210 Hz. There is 
also an offset chosen so that during the oscillation period 
it goes from zero to the maximum amplitude and then 
back to zero, always pushing and never pulling. The max- 
imum amplitude reached by the field gradient character- 
izes the perturbing strength and has been varied from 
zero to 190mG/cm along the vertical axis. 

Finally, the atoms are held trapped for another 20 ms, 
after the time-varying field is turned off. The measure- 
ments are performed via time-of-fiight absorption images, 
with the atomic cloud is expanding in free fall. For small 
amplitudes and short time lapses we observe dipolar, 
quadrupolar, and scissor modes but no vortices. By in- 
creasing the values of both parameters we were able to 
produce vortices, which grew in number with the ampli- 
tude and/or the excitation elapsed time. Further infor- 
mation can be found elsewhere P, [3, [3 113 ■ 

During the experimental runs, the clouds may un- 
dergo a few different excitation modes before the vor- 
tices are produced. To better observe the modes we kept 
the current in the coils oscillating around 170 Hz during 
texc = 100 ms. Then, we wait a time t lapse before re- 
leasing the sample from the trap, allowing it to freely 
expand for 23 ms. A resonant probe laser pulse of 70 /zs 
takes an absorption image. The RF forced evaporation 
was adjusted to keep the condensed fraction at about 
60 ms. The temperature of the cloud was kept around 
100 nK, for the conditions of this experiment. After the 
absorption image is acquired, we fit the images using a 
bimodal profile, constituted by a Thomas-Fermi and a 
Gaussian profiles for the condensed and the thermal frac- 
tions. Finally, different aspects of these two components 
are observed as a function of the hold time, r. 



III. OVERALL MOTION OF THE CLOUD 
DURING OSCILLATION 

To follow the cloud's motion into the trap, we have 
taken a sequence of TOF snapshots varying r. The su- 
perposed images were stacked together to better visual- 
ize the full path, as shown in Fig. [TJ The arrows point 
towards the increasing r, evolving clockwise. As the exci- 
tation progresses, a variety of collective modes take place 
absorbing the energy pumped into the cloud by the os- 




250Mm 



FIG. 1. Overview of the cloud's motion path after 23 ms of 
time-of-flight (TOF) absorption images. The arrows point 
towards the increasing hold time, t. 

dilating magnetic field. 

First, the clouds' center-of-mass oscillate inside the 
trap in the dipole mode pTj . This mode corresponds 
to the full oscillation, in trap, with the trap characteris- 
tic frequency in each direction. Besides the dipole mode, 
the long and short axes of the cloud also oscillate out- 
of-phase, characterizing the collective quadrupole mode 
of excitation . Since the longitudinal and transversal 
directions do not oscillate exactly 90° out-of-phase, we 
believe that higher order modes might also be present. 

Additionally to these collective modes, the tilting of 
the condensed cloud's long axis is also observed as pre- 
viously reported [l|- This tilting oscillation is related to 
scissors mode Q ^-nd shows that the excitation is able to 
produce rotation in the cloud. 

An interesting observation seen in Fig. [T] shows that 
the dipole mode trajectory does not close on itself which 
is an evidence that the dipole mode is also excited in 
the direction of the weak confinement (lower frequency). 
In fact, the cloud is under dipole excitation in all three 
directions, with different amplitudes, resulting in a more 
complex trajectory for this cyclic path. 

The excitation of the collective modes seems to be play 
an important role on nucleating vortices in our system 
[l[. By changing the excitation frequency we were able 
to vary the amplitude of oscillation of those described 
modes. We also observed a relative motion between the 
two components of the fluid. This out-of-phase oscil- 
lation in between the thermal cloud and the condensate 
correspond to a mode previously investigated in detail by 
Ketterle and van der Straten [13, [13 ■ It is quite surpris- 
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ing that this mode is coupled to the dipole mode and it 
may be fundamental to the vortices nucleation. The con- 
densed and the thermal fractions oscillate with respect 
to each other in a macroscopic scale such that there is 
a relative motion between the two centers of mass. Fig- 
ure [2] presents three arbitrary density profiles during the 
oscillation. The presence of such mode, excited by the 
oscillatory motion, is quite relevant to understand the 
mechanisms of vortex generation and evolution [l[ . 

The largest amplitude for the relative motion takes 
place along the cloud longest axis, but it is also present 
along the short axis since the excitation is not limited 
to a single axis. To determine the in situ sizes from the 
expanded cloud values, we calculated the scaling factors 
following the method developed in Ref. [1^. The results 
for the in situ counterflow motion are shown in Fig. [3] 
for the longitudinal and the transversal directions. The 
smallest displacement observed between the two clouds 
corresponds to the upper and lower turning points seen 
in Fig. [TJ In these points the dipole mode has the maxi- 
mum acceleration and the counterflow has the maximum 
velocity. 

While the longitudinal relative position in between the 
two clouds center of mass reaches distances of about 
30 /im, the transverse relative motion reaches the largest 
displacement of about 12 fim. Both motion occur in same 
frequency, as one can see on the plot in Fig. |3l 

We analyzed the motion of this out-of-phase mode us- 
ing the force per unit volume as derived by Gorter and 
Mellink [ll: 



(1) 



where and p„ are the densities of the superfluid and the 
normal fluid, respectively, vq is a constant on the order 
of 1 cm/s and A is a temperature dependent constant. 
It should be noted that Eq. ([T]) is a phenomenological 
relation, which is proportional to the cubic of the relative 
velocity (v„s = v„ — v^) in between the two fluids, and it 
seems to be valid even beyond the turbulence onset in the 
superfluid. The origin of this interaction is the scattering 
of thermal excitation as discussed by W. F. Vinen 2J|. 

The out-of-phase oscillatory motion corresponds to rel- 
ative velocities up to v„s w 3 cm/s under our experimen- 
tal conditions. This velocity can be larger or smaller, 
depending on the amplitude of the oscillating external 
magnetic field. 

To better understand the net result of the superfluid 
fraction slosh motion, we analyzed the surface ripples ap- 
pearing on the boundary region separating the condensed 
and thermal fraction (see Fig. [5]). The ripples may be de- 
scribed as undulations, or ruffling of the BEC/thermal in 
between surface giving it a wavy form and which may end 
up folded on itself like wave crests on the shore forming 
tubes or vortices. To analyze it, we have first fitted the 
thermal fraction only and then subtracted it from the 
image, resulting in a pure EEC shot, which was then 
fitted by the Thomas-Fermi profile and, finally, it also 




FIG. 2. Sequence of arbitrary picked density profiles showing 
the relative motion in between the thermal (normal fluid) and 
the condensed (superfluid) fractions. 



subtracted from the image, this resulted in a perimeter 
line surrounding the condensed fraction. 

To help on the data analysis we defined the surface rip- 
ples characteristic length. It is determined as the high- 
est peak resulting from the Fourier transform over the 
perimeter line separating the condensed and the thermal 
fractions on the sloshed clouds. In Fig. SI we show the 
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FIG. 3. Relative position between the centers of mass of the 
superfluid and thermal components for the longitudinal and 
transversal directions as a function of hold time after the ex- 
citation. The dotted lines are sinusoidal fittings. 
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FIG. 4. Characteristic length of the ripples in the superfluid 
surface as a function of the hold time. The dotted line is 
an eye guide only. For reference, we displayed the value of 
healing length. 



characteristic length of the ripples versus the hold time 
in order to allow for comparing to Fig. [S] We observe 
that the ripples amplitude are much larger around the 
turning points, where the maximum relative velocity in 
between the superfluid and normal fraction takes place. 

The larger the ripples the more likely to nucleate vor- 
tices into the superfluid fraction it will be during the 
slosh motion. To show these surface ripples, we present 
a typical time-of-flight cloud in Fig. [5] As it can be seen, 
the ripples may be the starting the vortex nucleation and 
one might be able to see the characteristic shape of the 
surface ripples already start to looking like vortices. In 
fact, a closer look reveals the vortex nucleation around 
the regions where the fluctuations are larger. 
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FIG. 5. Typical time-of-flight absorption image showing the 
ripples on the superfluid surface see the black perimeter region 
due to the counterflow. 



The experiment may be understood as a counterflow 
of the superfluid over the normal fraction. We believe 
that when the counterflow exceeds the critical velocity, 
the surface ripples grow larger, and may induce the for- 
mation of vortices. In our experiment due to the oscilla- 
tory nature of the counterflow, the maximum velocity 
is sustained for a short time lapse, with small vortex 
nucleation less likely evolving to the turbulent regime. 
Different experimental condition with larger amplitudes 
already demonstrated the onset of turbulence Q. 

Assuming the mutual friction force, Eq. ([1]), as the 
cause for the increasing fluctuations on the interface be- 
tween the two fractions, the regions with larger ripples 
correspond to Vns > vg. Looking at Fig. 21 for r running 
from 3.5 to 5.5 ms, v„s < vq and fluctuations are kept low 
(basically zero characteristic length). Around 5.5 ms, the 
relative velocity must be close to vq, and from there on 
the increase in Vns shall exceed vq. Hence, we estimate 
vq to be about 2.5cm/s. 

We have analyzed the situation using the arguments 
found in [15| , where the turbulence in superfluids is inves- 
tigated in terms of two dimensionless parameters: q and 
Rcs- Basically, the intrinsic parameter q deals with the 
relative importance of two competing terms determin- 
ing the energy dissipation in the system. The extrinsic 
parameter Rcg, the superfluid Reynolds number, which 
characterizes the flow of the superfluid fraction, given by 



Re, 



2TTh 



(2) 



where ^ is the healing length, v the flow velocity and m 
is the atomic mass. 

These parameters can be considered in the case of two 
mutually penetrating fluid components. Considering v = 
v„s and ^ = (87rpsas)~^/^, where a, is the s-wave scatter- 
ing length (~ lOOflo) and ps — 2.5 x 10^"' atoms/cm"^, we 
obtain Re, running from to 1. According to reference 
, these values are not large enough to generate turbu- 
lence. This statement is in agreement with our observa- 
tions, where the fluctuations are able to nucleate vortices 
but they do not evolve to a turbulent state. The normal 
component of the fluid vanishes close to T = K stopping 
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the counterflow and the frictional force. As a result, we 
should expect the vortex nucleation and the evolution of 
the vortex configuration to become more difficult as the 
temperature goes to zero. In fact, this justifies some of 
the previous observations [1^ as well as the need to in- 
troduce dissipation in the theoretical jpodels to be able 
to predict the turbulence emergence [l2|. 

IV. CONCLUSIONS 

We have investigated the slosh move of a non pure 
BEC sample, containing a condensed (superfluid) em- 
bedded in a thermal fraction (normal fluid), excited by 
time-varying magnetic fields. The relative motion be- 
tween the two components was interpreted as an out-of- 
phase oscillation mode, previously investigated by other 
groups [13, [13 • It was observed that ripples are more 



intense near the time instants of largest relative veloc- 
ity, where high amplitude ripples are generated on the 
superfluid/normal fluid interface, creating the right con- 
ditions for the vortex nucleation. Our findings helped 
on explaining previous experimental results [1], as well 
as the presented results may bring new possibilities for 
exploring the superfluid character of atomic BECs. 
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